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Разработана и реализована математическая модель нестационарных разделительных процессов обогащения урана в промы- 
шленных центробежных каскадах, которая может быть использована в компьютерном тренажере для подготовки специалистов 
разделительного производства и применения в качестве экспертной системы в автоматизированной системе управления техно- 
логической схемой. 


Введение 

Известно, что в общем случае нестационарные 
процессы в разделительном каскаде описываются 
нелинейным уравнением второго порядка в частных 
производных (относительно времени и номера сту- 
пени) параболического типа с переменными коэф- 
фициентами, аналитические решения которого по- 
строены лишь для ограниченного числа задач [1]. 
Значительно более перспективными методами ис- 
следования нестационарных процессов, являются 
методы численного решения уравнений описываю- 
щих эти процессы. Вопросам численного исследо- 
вания нестационарных процессов в разделительных 
каскадах посвящено много работ, например, [2]. Од- 
нако используемые в этих работах математические 
методы либо имеют ограниченную область практи- 
ческого применения, либо требуют больших вре- 
менных затрат на численное решение, либо не учи- 
тывают особенностей построения технологических 
схем центрифужного производства. При этом в ра- 
ботах рассматривается простейший симметричный 
противоточный прямоугольно-секционный каскад 
[3], в котором секции ступеней работают в одинако- 
вом гидравлическом режиме. Между тем, эффект 
обогащения центрифуги значительно выше (при- 
мерно в 100 раз), чем у диффузионных машин, а 
длина центрифужных каскадов существенно мень- 
ше. Как правило, ступени работают в отличающих- 
ся гидравлических режимах. В этих условиях повы- 
шается вклад каждой ступени в работу каскада, и со- 
ответственно при исследовании нестационарных 
процессов важно знать состояние каждой ступени. 

В работе предложена математическая модель, 
описывающая практически все виды нестационар- 
ных процессов как в отдельной ступени, так и в ка- 
скадах произвольного построения, отличающихся 
от традиционных симметричных противоточных ка- 
скадов. К тому же, моделируется вся многокаскад- 
ная технологическая схема, представляющая собой 
совокупность каскадов вместе со связывающими их 
линиями межкаскадных коммуникаций (МКК). 

Модель является дополнением к ранее разрабо- 
танной математической модели нестационарных ги- 
дравлических процессов в многокаскадной техноло- 
гической схеме (ТС) центрифужного завода [4, 5]. 


1. Описание математической модели нестационар- 
ных разделительных процессов 

Рассмотрена многокаскадная ТС разделитель- 
ного завода, состоящая из А- каскадов (рис. 1). 
Каждый каскад т (1 <т<К) состоит из ступеней с 
нумерацией от отвала к отбору. 

Схема соединения ступеней в каскаде может 
быть произвольной, например, так называемой па- 
раллельно-последовательной (рис. 2). 



Рис. 1. Многокаскадная ТС 


С„ 2 



Рис. 2. Параллельно-последовательная схема соединения 
ступеней в каскаде 

Ступень с номером п состоит из 5 т параллельно 
соединенных секций; секция і ш (\<і пт <5 ш ) состоит из 
•4,,, параллельно соединенных газовых центрифуг. Га- 
зовые центрифуги в секции считаются идентичны- 
ми, секции могут быть различными (разное число 
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ГЦ, разные геометрические размеры коммуникаций 
и т.д.). Конструктивные параметры К (порядка 
4.. .12), Ы т (порядка 4.. .25) определяются номенкла- 
турой выпускаемой продукцией, изменением ком- 
плектации и технологических характеристик ГЦ. 

В предложенной работе рассмотрена однообъе- 
мная модель разделительной ступени. Каждая сту- 
пень каскадов многокаскадной технологической 
схемы центрифужного производства считается за- 
полненной бинарной смесью с одинаковой концен- 
трацией изотопа 235 Ц, равной его концентрации в 
коллекторе питания этой ступени. Газосодержание 
разделительной ступени считается равным сумме 
газосодержаний всех объемов, входящих в состав 
ступени и учитываемых при расчете нестационар- 
ной гидравлики ТС [4]. На выходе из ступени (в кол- 
лекторах отбора и отвала) концентрации, определя- 
емые разделительными свойствами ступени, соот- 
ветственно выше и ниже концентрации питания. 

Основным уравнением нестационарного про- 
цесса разделения является уравнение баланса обо- 
гащаемого компонента питающей смеси в ступени 
п каскада т многокаскадной ТС: 

<Н м ^С 0т ) = Т т {і) + ^8 Яя Т Якя К Яя + 
сіі ы 

+ X ^пт р ^т р ^т р + 

к = 1 р=п + 1 \ѵ= 1 

+8 ОР С 3аК р„{0 Л Мт — ^ОЛи (О) — Х'^'Оіііюі^'Оіиі > (1) 

/=1 

где К - количество каскадов ТС; 5 - количество 
технологических секций в ступени я каскада т; 
М т - газосодержание ступени я каскада я?; С 0пт - 
концентрация питающей смеси ступени я каскада 
т; Т п „({), К 0п (Г) - расход и концентрация внешнего 
потока питания; 5 Пп - признак связи потока отбо- 
ра Т Шп с концентрацией К ш „ каскада к с питанием 
ступени я ( Ыт ) каскада т (к^т); 5 т „ - признак 
связи потока отвала Т п „ с концентрацией К Пп ка- 
скада к с питанием ступени я (кфт) каскада т 
( Ыт)\ 8„ Тр - признак связи потока отвала 0 Тр с кон- 
центрацией С Тр ступени я каскада т с питанием 
ступени я каскада я?; 8 пѢ - признак связи потока 
отбора с концентрацией С ѣ ступени \ѵ каскада 
т с питанием ступени я каскада /я; 8 0Р - признак 
отборной ступени (0 или 1); 0 Шт , Т 0Рт - соответ- 
ственно потоки отбора отборной ступени 7Ѵ„, 
( т=\...К) и каскада т; С 3трт - концентрация пото- 
ка закрутки каскада /я; Ь Шт - поток питания сек- 
ции / ступени я каскада т. 

В общем случае входящие в ур. (1) величины га- 
зосодержаний ступеней, потоков питания каскадов 
(внешних относительно ТС и межкаскадных пото- 
ков, являющихся потоками отбора и отвала других 
каскадов, а также межступенных потоков) являют- 
ся функциями времени. Эти функции могут быть 
получены решением гидродинамических уравне- 
ний, которые рассмотрены в работах [4, 5]. 


Концентрации межкаскадных потоков питания 
К ш и К Шп определяются через соответствующие 
концентрации потоков отбора и отвала каскадов с 
учетом транспортной задержки т, аі под которой по- 
нимается время прохождения гидравлического 
возмущения по соответствующим линиям МКК. 
Время т ы определяется из решения уравнений не- 
разрывности для установившегося одномерного 
потока, изолированного от внешней среды [6], от- 
носительно неизвестной скорости протекания газа 
V/ по /-ому участку линии МКК: 

С = 5,рр І , 

где С - расчетная величина потока в линии МКК 
(одинаковая для всех участков линии МКК), г/с; 
/?,■ - площадь сечения /-го участка линии МКК, см 2 ; 

цЯІЗЗЗ 

р, = — плотность газа в /-ом участке линии 

МКК, г/см 3 ; и - молярная масса газа, г/моль; 
Р; - расчетная величина давления газа в /-ом участке 
линии МКК, мм рт.ст.; Я, - универсальная газовая 
постоянная, эргДмольК); Т і - температура газа, К; 
V/ - скорость газа в /-ом участке линии МКК, см/сек. 

Величины давлений в і-х участках линии МКК 
определяются из принятого турбулентного харак- 
тера протекания газа, известного из решения ги- 
дравлической задачи значения какого-либо давле- 
ния (например, давления в коллекторе питания от- 
борной ступени) и заданных значений геометриче- 
ских размеров этих участков. 

Отборные ступени каскадов имеют некоторые 
особенности, связанные с тем, что часть потоков 
отбора через регуляторы закрутки поступает обрат- 
но в качестве питания этих ступеней. Кроме того, 
на части каскадов закруточный поток подается по- 
сле прохождения очистительного каскада (ОК). 
Следовательно, концентрация закруточного пото- 
ка С 3ащт изменяется с некоторой транспортной за- 
держкой по отношению к концентрации потока от- 
бора отборной ступени. Время транспортной за- 
держки также определяется в задаче нестационар- 
ной гидравлики [4, 5]. 

Ур. (1) составляется для каждой ступени каска- 
дов многокаскадной ТС. Полученная система ура- 
внений состоит из _ЮѴ т -дифференциальных уравне- 
ний с 3 КМ т неизвестными, каковыми являются 
концентрации ступенных потоков питания, отбора 
и отвала (С 0т , С Рт , С Тт ). Для понижения размерно- 
сти системы воспользуемся уравнением сохранения 
легкой компоненты и соотношением, определяю- 
щим полный коэффициент разделения ступени^,,, 

г =0 С„ +(0 -1)0. , (2) 

О пт пт Рпт 'пт ' Тпт 9 \ / 


где Ѳ т - коэффициент деления потоков ( Ѳ ІШ = ) . 


О п . 


Хит 


С Ртп О С Тпт ) 

С Тпт (1 “ С Рпт ) 


( 3 ) 
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Из ур. (2) и (3) можно определить С Тш : 


г 

^ Тпт 


с„ 


X пт Ср 1т (х„ т 1) 


(4) 


Путем подстановки выражений для С Тт и С 0ш 
система ур. ( 1 ) сводится к системе дифференциаль- 
ных уравнений относительно Ш т неизвестных С Рт 
как функций времени. Решение этой системы ура- 
внений при заданных начальных значениях С Рт 
моделирует нестационарные разделительные про- 
цессы в каскадах ТС. 


2. Алгоритм решения уравнений нестационарного 
разделения 

Обозначив правую часть ур. ( 1) после подстанов- 
ки в него выражений (2) и (4) через у/ т (1) получим: 

^1 4Т„,„ (Ѳ ІШ ,С Рпт - ( 1 - Ѳ тп )/ Сі (С Рпт ) ) ] _ ^ ^ 

ііі 


где /с,( С Рш ) - функция от концентрации отбора сту- 
пени п каскада т , определяемая выражением (4). 

Система уравнений (5) представляет собой за- 
дачу Коши с заданными начальными условиями. 
Задача Коши за редким исключением не имеет 
аналитических методов решения. В численном ме- 
тоде вместо отыскания непрерывных функций вре- 
мени рассчитываются значения искомых функций 
в дискретные моменты времени, а система диффе- 
ренциальных уравнений заменяется тем или иным 
способом разностными уравнениями. В применяе- 
мом алгоритме, аналогично гидравлической задаче 
[4, 5], используется неявная схема Эйлера с перес- 
четом [7], которая обладает третьим порядком точ- 
ности на шаге и вторым на интервале. В соответ- 
ствии с ней, имеем: 

А/Г (К+1) т (к+1) г (к+1) _/і _д(*+1)\ Лк+\) (Мк+ 1) \\ _ 

1Ѵ1 пт \°пт ^ Рпт V 1 °пт /АО УУРпт )) 

м'Жс™ -а -О/АС,! )) = 

= уЛ +І, +'0- (6) 

Написав уравнения (6) для всех ступеней каска- 
дов ТС, получим систему алгебраических нелиней- 
ных уравнений относительно С^ 1 ’. Нелинейность 
обусловлена наличием членов, содержащих выраже- 
ние На каждом временном шаге проводит- 

ся линеаризация её нелинейных членов. В предло- 
женном алгоритме она заключается в замене функ- 
ции /сі +,І С ( Рт ’ уравнением касательной в точке С Шт . 

С ( Т К+1) =а^ +Ы Т к) С 1* +1) . (7) 

Тпт Тпт Тпт Рпт V / 


Коэффициенты а$ т и Щ ] пт вычисляются из условия 
равенства линейных (7) и нелинейных (4) выражений, 
а также их производных в некоторой точке С Шт . 


г 

^ОРпт 

X пт ^ 0 Рпт (X пт 


а Тпт + ^ТппТ'йРпт ’ 


( 8 ) 


(0}іт) _ 


\-Х п/и ^ОРпт^Х пт ^ 1 Г 


= Ьп 


(9) 


Из решения ур. (8, 9) получаем: 

Тпт \ у*- К) - С ( ~ К) ( у І ~ К) -П1 2 ’ 

I -А, пт ^ Рпт пт 

У (К) 

и(к) А , пт 

Т"т ~ \ -С ІК) (У ІК) -V)? ' 

После линеаризации получим систему линейных 
уравнений относительно неизвестных О р *'\ являю- 
щихся приближением к величинам С^ІІ+АІ). Си- 
стема решается методом исключения Гаусса [7, 8]. В 
качестве начальных условий, исходные значения 
Сй я =С Ди (И)), принимаются расчетные данные ис- 
ходного стационарного режима исследуемой ТС. 

Значения коэффициентов полного разделения 
Хт ступеней определяются разделительными свой- 
ствами ступеней. Разделительные свойства отдель- 
ной центрифуги (ГЦ) описываются эмпирической 
зависимостью, определяемой в результате матема- 
тической обработки экспериментальных данных, 
полученных при испытаниях оборудования данно- 
го типа ГЦ. 

х пт = ХппМ,Ѳ), (Ю) 

где §, Ѳ - расход питания и коэффициент деления 
потоков отдельной ГЦ. 

Считая величину х инвариантом, независящим 
от концентрации потоков питания ступеней, мож- 
но связать величины С Тт и С Рш по формуле (4) при 
работе ступеней в области концентраций, отли- 
чающихся от концентрации природного сырья. 

Разделительная характеристика Хт(8>Ѳ) опреде- 
ляется при работе ГЦ в стационарном режиме, ког- 
да выполняется баланс вещества и легкой компо- 
ненты. Использование этой характеристики в не- 
стационарном процессе объясняется тем, что про- 
цессы установления равновесия в отдельной ГЦ 
протекают значительно быстрее, чем в каскадах 
ТС. Центрифуга рассматривается как работающая 
в стационарном режиме, который медленно изме- 
няется. 

По описанной модели в среде Вогіапсі Оеіріті 5.0 
реализован программный комплекс и проведены 
расчеты для следующих случаев нестационарных 
возмущений ТС: изменение потоков отбора и пита- 
ния; перенос точек подачи питания; закрытие сту- 
пеней и секций каскадов; изменение частоты пи- 
тающего тока ГЦ; изменение загрузки оборудова- 
ния; изменение концентрации внешних потоков 
питания. В качестве разделительной использова- 
лась приближенная модельная характеристика, от- 
личающаяся от фактической, но отражающая ос- 
новные тенденции протекания реального процесса. 

В качестве примера, приведем результаты чи- 
сленного расчета нестационарных процессов раз- 
деления, возникающих при закрытии отдельных 
ступеней девятиступенного каскада (рис. 3). 

Сравнительный анализ пар кривых 1, 2 и 3, 4 
позволяет сделать вывод, что созданная модель 
адекватно описывает поведение реального объекта. 
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Повышение со временем обогащения отбора 
каскада (1,2) объясняется тем, что с момента зак- 
рытия ступени каскад переходит в гидравлический 
режим с повышенной нагрузкой на оставшиеся 
ступени, при котором величина % имеет большее 
значение. Ломаный характер кривых 2, 4 объясня- 
ется способом измерения концентрации 235 П, осно- 
ванным на регистрации у-излучения рабочей газо- 
вой смеси, интенсивность которого носит вероят- 
ностный характер. 





Рис. 3 . Временные зависимости отклонений концентрации 
отбора каскада относительно исходной в условных 
единицах (1 у.е. =0,001 %). Закрытие ступеней: 1) рас- 
чет; 2) факт ~ отборной и 3) расчет; 4) факт - предо- 
твальной 

Из рисунка также следует, что влияние отбор- 
ной ступени на конечный продукт значительно вы- 
ше влияния отвальных ступеней. 

Проведенные расчеты показали, что характери- 
стики быстродействия модели позволяют опера- 
тивно производить расчеты многокаскадной тех- 
нологической схемы в режиме реального времени. 
Так, например, время счета 4-х часового временно- 
го интервала протекания нестационарного процес- 
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са для четырехкаскадной ТС с общим количеством 
ступеней порядка 50 и временным шагом, равным 
1 с, составляет примерно 140 с. В связи с этим раз- 
работанная модель нашла применение в АСУТС в 
качестве экспертной системы. Кроме того, разра- 
ботанная модель используется в компьютерном 
тренажере для подготовки специалистов раздели- 
тельного производства [9, 10]. 

Заключение 

Предложена математическая модель нестацио- 
нарных разделительных процессов в промышлен- 
ных центрифужных каскадах, учитывающая отли- 
чительные особенности построения исследуемой 
технологической схемы разделительных произ- 
водств, которая представляет собой систему диф- 
ференциальных уравнений первого порядка с за- 
данными начальными условиями. Разработан алго- 
ритм решения уравнений, описывающих неста- 
ционарные процессы в центрифужных каскадах. 

Реализованная модель позволяет получать ре- 
шения в режиме реального времени при произ- 
вольном наборе возмущающих воздействий и ис- 
следовать нестационарные процессы, возника- 
ющие при изменении величин внешних потоков, 
их отключений; изменении частоты и отключении 
питающего тока; изменении точки подачи пита- 
ния; отказах и перенастройке межступенных регу- 
ляторов; отключении части оборудования; наруше- 
ниях в работе оборудования линий МКК; при раз- 
личных сочетаниях этих воздействий. 

Модель внедрена и используется в АСУТС раз- 
делительного производства ФГУП ПО "Электрохи- 
мический завод" в качестве экспертной системы, а 
также в компьютерном тренажере для подготовки 
специалистов разделительного производства. 
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